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Solar cells are engines converting energy supplied by the photon flux into work. All known 
types of macroscopic engines and turbines are also self-oscillating systems which yield a periodic 
motion at the expense of a usually non-periodic source of energy. The very dehnition of work in 
the formalism of quantum open systems suggests the hypothesis that the oscillating “piston” is a 
necessary ingredient of the work extraction process. This aspect of solar cell operation is absent in 
the existing descriptions and the main goal of this paper is to show that plasma oscillations provide 
the physical implementation of a piston. 


INTRODUCTION 

The standard model of work extraction is based on an engine composed of a working medium, a piston and two heat 
baths that are at equilibrium with different temperatures. Its importance stems from its success to set an universal 
bound to any work extraction process, the Carnot bound, which shows, in agreement with the Kelvin’s formulation 
of the second law of thermodynamics, two different temperatures are needed for work extraction. Besides, there is a 
complementary picture of an engine as a self-oscillating system “focusing on their ability to convert energy inputted at 
one freguency (usually zero) into work outputted at another, well-defined frequency” [1]. This picture seems to apply 
to all types of turbines and motors [1, 2] and one can expect that it is equally correct for engines powered by a flux 
of photons. 

The standard description of the photovoltaic cell involves the following processes [3, 4]: 

i) generation of the charge carriers due to the absorption of photons, 

ii) subsequent separation of the photo-generated charge carriers in the junction. 

As noticed in [4] the often used explanation of the second process as caused by the emerging electric field in p-n 
junction cannot be correct. Charge separation is supposed to produce a DC current which, on the other hand, cannot 
be driven in a closed circuit by a purely electrical potential difference. A standard thermodynamical explanation of 
electric current (work) generation in photovoltaic and thermoelectric heat devices is the following : 

Electrons gain energy in a form of heat current Jh from the hot bath, then flow against potential difference <I> 
producing useful power P = Je^, where Je is an electric current. A part of the heat described by the heat current 
Jc is dumped to the cold bath. 

The laws of thermodynamics put the following constraints 

Je^ = Jh — Jc (1) 


Jh 



( 2 ) 


Adding kinetic equations describing the processes of electron-hole creation, thermalization to the ambient temperature 
and recombination one obtains the correct formulas for the open circuit voltage and voltage-current relation. However, 
this picture does not explain the mechanism of persistent steady work extraction. Similarly, for a steam engine the 
net pressure due to the temperature difference obviously provides the net force acting on the piston but to explain 
the permanent periodic action of this engine we have to understand the details of operation of a piston linked to a 
flywheel and valves. 

This is exactly the place where the mechanism of self-oscillation supported by the external constant energy flow 
enters into the game. In the following we propose a model in which plasma oscillations play a role of the periodic 
motion of a “piston” and show that, indeed, under realistic assumptions a positive power is supplied to this essentially 
classical macroscopic oscillator. Subsequently, the collective charge oscillations at THz frequencies are rectified by a 
p-n junction diode (“valve”) to the output DC current. 

The mathematical formalism is based on the quantum Markovian master equations for slowly driven open quantum 
systems studied in [5] (compare [6], [7] for the fast driving case), and consistent with thermodynamics. 
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MODEL OF QUANTUM ENGINE 

We consider a model of heat engine which consists of a “working medium” called simply a system, two baths at 
different temperatures, and a “work reservoir” called often a “piston” which is a system supplying to or extracting 
work from the working medium. Because work, in contrast to heat, is an ordered and deterministic form of energy 
we expect that a piston should be a macroscopic system operating in the semi-classical regime. Therefore, within the 
reasonable approximation can be replaced by external deterministic driving [8] . 

For the readers convenience we briefly present the formalism of Quantum Master Equations (QME) based on the 
Davies weak coupling limit [11], the Lindblad-Gorini-Kossakowski-Sudarshan generators and their extension to slowly 
varying external driving. Then the thermodynamical consequences are discussed and a special generic class of models 
with diagonal, weak driving, applicable to the theory of solar cells, is presented. 


Master equation for open system with constant Hamiltonian 


The total system consisting of a system S weakly coupled to a bath B. The total Hamiltonian is a sum of three 
terms: 

i) a ’’bare” Hamiltonian of the system - replaced in the final formulas hy H - a“ physical renormalized Hamiltonian” 
containing the lowest order corrections, 

ii) Hamiltonian of the bath - iJg, 

iii) a system-bath (elementary) coupling: 

H,^t = A®F , {F)b = 0. (3) 

where A and F are observables of the system and the bath, respectively. (E)g denotes the average with respect to 
the stationary state of the bath. 

Two main ingredients enter the QME derived using Davies weak coupling limit procedure [11]: 
a) the spectral density of the bath 

/ + 00 

e—‘(F(t)F)Bdt>0, (4) 

-OO 

b) Fourier components of the coupling operator 

A{t) = U\t)AU{t) = U{t) = (5) 

{^} 

Introducing system Hamiltonian spectral decomposition and Bohr frequencies 


one obtains the relations 


H = 'Y^k\k){k, {iv} = {{tk - ei)/h} 

k 


( 6 ) 


[F[,A{uj)] = —hujA{ui) , A{—uj) = A\uj). (7) 

where A{ijj) called transition operators or Lindblad operators correspond to energy exchange of hw. The standard 
derivation yields the QME in the Schrodinger picture 

^ = -^-[H,pit)]+Cpit) (8) 

where 

£p= i^G(^)([A(a.)p,At(a.)] + [A{u;), pAHu;)]) (9) 

Eor general interactions Hint = ® 

Cp = i^^Ga/3(w)([A„(w)p,A^(w)] -f [A„(w),pA^(a;)]) 

a,/3 {i^} 

where [Gapioi)] is a positively defined relaxation matrix. 


( 10 ) 
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Properties of QME 

The QME’s obtained in the weak coupling limit possess the following properties [12], [13]: 

1) £ possesses the Lindblad-Gorini-Kossakowski-Sudarshan structure what implies that the solution of QME is a 
completely positive, trace preserving, one-parameter semigroup. 

2) Hamiltonian part commutes with the dissipative one, i.e. 

p{t) =U{t)e*^p{Q) = e*^U{t)p{Q), U{t)p = U{t)pU\t). (11) 

3) Diagonal (in H - basis) and off-diagonal density matrix elements evolve independently. 

4) The stationary state p satisfying [H, p] = 0 and Cp = 0, always exists for finite-dimensional Hilbert spaces. 

5) For a single heat bath (i.e. a reservoir in the thermal equilibrium state): 
a) spectral density satisfies the KMS condition 


b) the Gibbs state is stationary 

G{-uj) = e-^‘^/'^^'^G{oj), 

(12) 


p^^ = Z-iexp{-^}, £p'=’=0. 

(13) 

c) the following quadratic form 

{X,C*Y)eq = Ti{p^'iX^C*Y) 

is negatively defined [12],[13] , where C* is the Heisenberg picture generator given by 

(14) 

C*X = 

-^^Ga/3(w)(A^[Ar, Aa(a;)] -1- ([A^,X]Aa(w)). 

(15) 


a,0 {w} 


It means that the Heisenberg picture generator C* can be treated as a hermitian, negatively defined operator acting 
on the space of (complex) observables equipped with the scalar product {X,Y)eq = Tr(p®«Xly). 


Entropy balance and the Laws of Thermodynamics 

We identify the physical entropy with the von Neumann entropy of the reduced density matrix S{p) = —ksTYip In p) 
and use also the relative entropy S{p\a) = A:BTr(plnp — plncr). 

For the solution p(t) of MME, and the stationary state p, {S{t) = S{p{t))) 

^S{t) = K{t) - kB-^Tr{p{t) Inp) (16) 

where K{t) = —-^S{p{t)\p) = —/cBTr([£p(t)][lnp(t) — Inp]) > 0 is interpreted as an entropy production. Positivity of 
the entropy production follows from the fact that for any completely positive and trace-preserving map A, S'(ApjAcr) < 
S'(pjcr) [14]. For many independent heat baths one obtains the Second Law in the following form 


dt 


y->o, 

VTfc - ’ 


where Jfe is a heat current flowing from the fc-th bath. 


(17) 


The case of a slow piston 

In order to define work we introduce the time-dependent Hamiltonian H (t) = Hq V{t) with slowly varying 
and typically periodic perturbation V (t) which gives a semi-classical description of a piston. We combine now the 
weak coupling assumption concerning the interaction of the system with several baths with a kind of adiabatic 
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approximation concerning the time-dependent driving. The former condition means, practically, that the relaxation 
rates are much smaller than the corresponding Bohr frequencies. The later is valid for the case when the time scale of 
driving is much slower than the time scale determined by the relevant Bohr frequencies. This is a similar situation to 
standard adiabatic theorem in quantum mechanics and implies that in the derivation of Master equation we can put 
the temporal values of Bohr frequencies {a;(t)} and transition operators {Aa(a;(t))} satisfying (7) with H replaced 
by H{t). 

Under the conditions of above one obtains the following form of QME 

(18) 

3 

where Cj{t) is a LGKS generator obtained by a weak coupling to the j-th bath and for a fixed H{t). The properties 
of generators Cj{t) imply the Zero-th Law of Thermodynamics {Pj = l/ksTj) 

e-PjHit) 

= = (19) 

Using the definitions [5]: lU-work provided by 5, Q “ heat absorbed hy S, E - internal energy of S 

E{t) = Tr(p(t)i7(t)) (20) 

- («) 

|<3(,)=T,(^H(i)) (22, 

= ^Tr(i7(t)£,(t)p(t)) = ^ ^Q,(t), (23) 

3 3 

where Qj is a heat absorbed by S from Bj^ one obtains the First Law of Thermodynamics 

The Second Law of Thermodynamics follows again from Spohn innequality 

fs{t) - Y, W ^ 0 (25) 

3 ^ 3 

where <jj (t) is an entropy production caused by Bj and given by 

aj{t) = kBTr{Cj{t)p{t)[\np{t) - lnp®'^(t)]) > 0. (26) 


Weak, diagonal and periodic driving 

We consider a generic case of oscillating weak driving V(t) which in the lowest order approximation can be replaced 
by the diagonal (in the basis of Hq) operator 

V{t) = gM sin Ut, [Hq, M] = 0, (27) 

with the small coupling constant g «1. In this case all Hamiltonians H{t) commute. 

The unitary part of the dynamics U{t) governed by H{t) commutes with Hq and M. One can write C{t) = £[5(<)] 
where £[^] is computed with the system Hamiltonian Hq -|- and ^(t) = gsuiLlt. Again for different ^’s the super¬ 
operators £[^] commute with the Hamiltonian part —i[H{t),-] and one can use their lowest order expansion with 
respect to ^ 


=^[0]+^£'[0] + O(e') 


(28) 
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In the next step we use the lowest order expression for the dissipative part of the super-propagator 

Aoit) = exp{ J £[^(s)]ds} 

~ +g [ (sinOs) ds. 

Jo 


(29) 


Applying now the definition of work (21), using the commutation properties of the generator and the fact that 
[U{t),M] = 0, one obtains the formula for the stationary average power output 


to->00 tg Jg at 


dt 


(30) 




= — 5 O lim — / 
to^oo to Jo 


Ti^MAn{t)p{0)^ cosOtdt. 


Inserting the lowest order expression (29) into (30) and using the fact that all super-operators commute one can 
compute the super-operator-valued integral like a usual one. Then, we apply the obvious assumption that the unper¬ 
turbed dynamics drives, asymptotically, any initial state p(0) to the stationary state denoted by p[0]. Finally, 
the limit tg — ?► 00 in (30) can be performed leading to the second order approximation for the average output power 


P = -g^Tifu 


02 


02 + (£[ 0 ])' 




'[0]p[0]). 


£[^] also possesses the stationary state p[^], i.e. £[^]p[^] = 0 and hence we can use the identity 

c'lm] = -mp'm 


(31) 


(32) 


where “prime” denotes the derivative with respect to Then, replacing the Schroedinger picture generator £[0] by 
the Heisenberg picture one £*[0] we can transform the formula (31) into 

P = -ls’Tr(p'[01jp^A_£.|0]M). (33) 

For the case when the decay rate of M is much lower than the modulation frequency we can neglect (£*[0])2 in (33) 
to obtain the simplified expression 

^ = -^52Tr(p'[0]£*[0]M). (34) 

The compact formula of above will be used to derive the specific expression for the solar cell power. 


No output power from a single heat bath 


The obtained lowest order formulas for power (33), (34) are still consistent with thermodynamics. Namely, assuming 
that the reservoir is a thermal equilibrium bath at the temperature T we have the following properties: 

1) p[^] is the Gibbs state with respect to the Hamiltonian Hq + ^M, 


2) the Heisenberg picture generator £*[0] is a negatively defined operator on the Hilbert space equipped with the 
scalar product {X,Y)eq = Tr(p[0]ArlF) (compare with (14), (15)). 

Using the first property one can compute p^[0] and rewrite the formula (33) as 

5^ 


2kBT 


2kBT 



(35) 

H2 + (£*[0])2'^ [0]M)e, <0 

(36) 


which is obviously negative as well as its simplified version (34). Therefore, as expected, one cannot extract power 
from a single heat bath. 
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Feed-back mechanism 


For a reservoir composed of two equilibrium ones at different temperatures a positive output power P > 0 can 
be obtained and in this case the mechanism of self-oscillations works. Treating the external perturbation as caused 
by a coupling of the system to a macroscopic oscillator we see that the positive power is supplied to the oscillator 
increasing the amplitude of its oscillation (positive feed-back) up to the moment when the net energy supply from 
the reservoir is compensated by the load attached to the oscillator. 


MODEL OF SEMICONDUCTOR SOLAR CELL 


A solar cell is an engine which produces work from heat exchanged with a non-equilibrium bath. The bath consists 
of the photonic non-equilibrium reservoir characterized by the local state population n[uj\ and the basically phononic 
heat bath at the temperature T of the device. The typical semiconductor solar cell consists of a moderately doped 
p-type absorber, on both sides of which a highly doped layer is formed, n-type on the top side and p-type on the back 
side, respectively. The electronic states in the valence band and in the conduction band are labeled by the index k 
which corresponds to the quasi-momentum /ik (spin can be easily added) with the energies £'«(k) and £’c(k) ■ We 
assume a direct band structure with vertical optical transitions which preserve quasi-momentum (see fig. 1). 



conduction 

band 


valence 

band 




k 


FIG. 1. Schematic picture of leading processes involving electrons, photons and phonons in a semiconductor with a direct band. 


The basic irreversible processes are the following: 

i) fast intraband thermalization processes mediated by phonons and described by Cth- 

ii) optical transitions between the valence and conduction band, which create or annihilate electron-hole pairs, 
described by Cem, 

iii) non-radiative electron-hole recombination which is neglected in our idealized model. 
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Plasma oscillations and current rectification 

The fundamental question in the presented approach to work generation in solar cell is the origin of periodic 
oscillations which can be seen as classical. The frequency 17 is assumed in our derivations to be much smaller than 
the frequency wq, but much larger than the recombination rate in order to justify (34). The only phenomenon which 
satisfies all these requirements is plasma oscillation visible for p-n junctions. Their appearance is due to the fact 
that a p-n junction creates an interface between regions of different electron concentrations which can oscillate in 
space producing collective macroscopic electric field oscillations. In several experiments such oscillations have been 
observed [9, 10], with typical frequencies 17/27r ~ ITHz, much lower than wq — IdOOTiJz corresponding to the energy 
gap ~ leV. On the other hand 17 is much higher than the recombination rate ^ what justifies the transition 

from (33) to (34). In the final step of cell operation the THz plasma oscillations must be converted into a direct 



FIG. 2. Schematic picture of rectification of plasma oscillations. Plasma oscillations caused by the collective motion of 
free carriers represented by the dimensionless amplitude ^(t) in the eq.(38). For small ^ plasma self-oscillation is described by 
the quantum harmonic oscillator coherently pumped by the feed-back mechanism. Due to the asymmetric diode-type potential 
at the edges, charge oscillations are converted into a DC current. 

current. A qualitative picture of this mechanism is shown on Fig. 2. The plasmonic degree of freedom is represented 
by the quantum levels in the asymmetric potential which is harmonic for lower energies. Asymmetry is due to the 
p-n junction which defines an “easy” direction for the carrier flow (to the left). The work supplied to the oscillator 
drives the unidirectional electric current. 


Hamiltonians, Master equations and stationary states 

The electrons in a semiconductor occupying the conduction and valence bands are described by the annihilation and 
creation operators Ck, cj^. and Wk: ^^k’ respectively, subject to canonical anticommutation relations. The unperturbed 
Hamiltonian reads 

Hq = ^(^Ec(k)cj^Ck + (37) 

k 

In the p-n junction a non-homogeneous free carrier distribution created in a self-consistent build-in potential can be 
perturbed producing collective plasma oscillations with the frequency 17. The associated time-dependent perturbation 
added to the electronic Hamiltonian (37) has a mean-field form [N - number of atoms in the sample) 

^(7)M = ^(7)-^F;g^(ctck(38) 

where ^ is a small dimensionless parameter describing the magnitude of deformation, Eg is the relevant energy scale, 
and cj^Ck, 'Ck'i’k are number operators of free electrons and holes, respectively. 
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To apply the formulas derived in the previous sections we notice first that the driving perturbation (38) depends 
only on the total numbers of both types of carriers and hence does not interfere with intraband transitions. Therefore, 
the relevant Bohr frequency is associated with the gap Eg yielding the time scale ~ much faster than the 

modulation period ~ This justifies the adiabatic approximation. The weak coupling assumption is for sure 

satisfied for very slow radiation recombination processes. 

Among the basic irreversible processes the intraband thermal relaxation is the fastest (thermalization time ~ 
and therefore, the stationary state of the electronic systems with the total Hamiltonian Hq+^M is, within a reasonable 
approximation, a product of grand canonical ensembles for electrons in conduction and valence band with the same 
temperature T of the device and different electro-chemical potentials p,c and respectively. The associated density 
matrix has form 



The electro-chemical potentials are determined by the numbers of carriers and hence by doping and radiative and 
non-radiative processes of electron-hole creation and recombination. 

Because the intraband thermalization to the ambient temperature T does not change the number of free electrons 
and holes, i.e. Ch^M = 0, the generator does not enter the formula for power (34). Here, one can doubt whether 
for such fast relaxation the weak coupling condition and hence the validity of the Markovian approximation leading to 
Cth holds. However, intraband relaxation does not contribute to work generation but only determines the structure 
of the stationary state. The form of this state expressed in terms of Fermi-Dirac distributions is generally accepted in 
the literature [4] and the accuracy of the Markovian approximation for the thermalization process is not very relevant. 
Finally, the contribution which remains in the eq. (34) describes the quasi-momentum preserving (vertical) transitions 
and reads 


CemPiP 



k 


+ 7ex(k)([cj^Uk, + [Ck^k P, 4ck 



(40) 


7T‘ec(k) 


— [l -I- n(a;k)l, 7ex(k) = —^(wk) 

Tse Tse 


(41) 


where Tse is the spontaneous emission time, fiwk = Ec{^) — Ey{h), and n{uj) denotes a number of photons occupying 
a state with the frequency w. 


Power and efficiency 

One can insert all elements computed in the previous section into the expression for power (34). Then we use 
the properties of the quasi-free (fermionic Gaussian) stationary state (39) which allow to reduce the averages of 
even products of annihilation and creation fermionic operators into sums of products of the Fermi-Dirac distribution 
functions 


/c(k) 


1 

gP(£;c(k)-/ic) _|_ 1 ’ 


fv{k) 


1 

g/3(£;„(k)-/xG + 1’ 


(42) 


with /3 = l/ksT. 

The leading order contribution to power possesses a following form 

k 

- 7rec(k)[l - /„(k)]/c(k)j, 


(43) 
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where Near = X)k [(‘^k''k)o + (vk?^k)o] is the total number of free charge carriers. 
One can introduce the local temperature of light T[w] defined by 

1 + n{Lo) 


For the incident sunlight on Earth a rough approximation holds 


l^sun (^) 


A 

^huj/ksTs _ \ ’ 


(44) 


(45) 


where ~ 6000/F is the temperature of the Sun surface and A = [Rsun/RoY — 2 x 10“® is the geometrical factor {Rsun 
- Sun radius) which takes into account the photon density reduction at large distance from the source. In particular, 
for the typical value of the energy gap Eg = hujo ~ leV the effective temperature of sunlight T[a;o] — lOOOiF. 

Because the product of population numbers for free carriers given by [l — /„(k)]/c(k) is essentially concentrated 
on the interval Wk € [wojWq + ©(fc^T/^i)], and T ~ 300iF <C T[ujo] ~ lOOOiF <C huJo/kg ~ 12000K, the expression 
(43) can be approximated by 


P = 


g'^ElNearF 
ksT Tse 



T 1 



(46) 


where F = + n(a;k)] [l — /j,(k)]/c(k) > 0 is independent of the size of a cell, and the voltage $ is identified 

with the difference of electro-chemical potentials, i.e. e$ = /ic — 

The condition for work generation by the solar cell reads 


ed) < e<I>o = rjc Eg, 


?7C = 1 - 


T 


T[u;o] 


(47) 


what implies that $o is an open-circuit voltage of the cell for the idealized case [15]. 

The presence of the Carnot factor gc suggests also the interpretation of the eq. (47) in terms of thermodynamical 
efficiency. Indeed, the incident photon of the frequency uj > ujq produces an excitation of the energy close to Eg in the 
process of electron-hole creation followed by the fast thermalization of an electron to the bottom of the conduction 
band, and a hole to the top of the valence one. Then, a part of energy Eg is transformed into useful work, equal at 
most e$o per single electron flowing in the external circuit. The maximal efficiency r]rnax under the conditions that 
each photon with the energy higher than the gap produces an electron-hole pair and non-radiative recombination 
processes are neglected, is given by the product rjrnax = Vu'Vc^ where rju is the so-called ultimate efficiency computed 
under the assumptions: 

a) “... photons with energy greater than Eg produce precisely the same effect as photons of energy Eg, while photons 
of lower energy will produce no effect” [17], 

b) the whole Eg is transformed into work. 

Under standard illumination conditions the ultimate efficiency of a solar cell can reach 44% and the Carnot factor is 
about 70% what yields gmax — 31% - the Shockley’s detailed balance limit [17]. Actually, photons are absorbed along 
their path in the absorber and n(uJo) decays exponentially with the penetration distance. Taking a more realistic 
average value T = (r[a;o] -I- T)/2 ~ 650K one obtains Pmax — 24%, much closer to the real performance of standard 
GaAs solar cells. 


Conclusions 

The presented model based on the idea of self-oscillations explains the dynamical origin of work generation in 
photovoltaic cells which is not present in the standard “static” picture. The main new ingredient is the role of plasma 
oscillation as a “piston” which transforms the steady heat input from the photon flux into periodic motion. This 
model provides a bridge between the theory of driven quantum open systems applied to heat engines and the theory 
of photovoltaic devices. The formulas (44) and (47) explain in a simple way the meaning of the “light temperature”, 
the Carnot bound, and the linear relation between the open circuit voltage and the device temperature. 

The experimental verification of this model should provide the evidence of THz plasma oscillation in the device 
with the amplitude square proportional to the power output. Such oscillations produce a weak THz radiation which, 
in principle, could be detected. 
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The similar ideas can be applied to other types of heat engines with “hidden self-oscillations”. It seems that 
thermoelectric devices based either on bimetallic or semiconductor p-n junctions can be described by the very similar 
models. Plasma oscillation remains a piston and sunlight is replaced by the hot bath. For organic photovoltaic 
systems, proton pumps or photosynthesis there exists quite strong evidence of the important role of coherent molecular 
oscillations played in the energy and charge transfer (see e.g. [18]). It is plausible that those oscillations can play the 
role of a piston in the work extraction mechanism as well. 
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